Nucleon-Nucleon Scattering Observables in Large-iV c QCD 
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Nucleon-nucleon scattering observables are considered in the context of the large N c limit of 
QCD for initial states with moderately high momenta (p ~ N c ). The scattering is studied in the 
framework of the time-dependent mean-field approximation. We focus on the dependence of those 
observables on the spin and isospin of the initial state which may be computed using time-dependent 
mean-field theory. We show that, up to corrections, all such observables must be invariant under 
simultaneous spin and isospin flips (i.e. rotations through n/2 in both spin and isospin) acting on 
either particle. All observables of this class obtained from spin unpolarized measurements must be 
isospin independent up to 1/N C corrections. Moreover, it can be shown that the leading correction 
is of relative order 1/N% rather than 1/N C . 



The strong interaction between two nucleons is the 
basic ingredient of nuclear physics. We wish to explore 
qualitative features of the nucleon-nucleon interaction 
which may be understood from QCD. In this context, 
it may be useful to consider the interaction in the limit 
when the number of colors, N c , of QCD becomes large 
and to treat l/N c as an expansion parameter. Cer- 
tain aspects of QCD can be deduced in this limit in a 
model-independent way. For example, the spin-flavor 
structure of certain amplitudes in the single baryon sec- 
tor may be fixed S. The large N c limit can also be 
used to determine the leading spin-isospin dependence 
of certain nucleon-nucleon scattering observables. 

Baryons in large- N c QCD were first discussed by Wit- 
ten, 0] , who argued that they are well approximated by 
a Hartree-type mean-field treatment. This picture works 
cleanly when dynamical gluons are integrated out. For 
the case of light quarks, explicitly deducing the Hartree 
equations of motion is not a tractable problem with 
present techniques. However, the N c dependence of var- 
ious baryon properties can still be deduced. Thus, for 
example, the baryon mass is of order N c and the baryon 
size is of order unity . 

Witten also considered the baryon-baryon interaction. 
He argued that the strength of this interaction is of or- 
der N c and the proper framework is the time-dependent 
mean- field theory (TDMFT). The argument is similar to 
the single baryon case: each quark moves in an average 
potential due to the other 2 N c — 1 quarks which in this 
context is time-dependent. The mean- field nature of the 
nucleon-nucleon scattering in the large N c limit imposes 
constraints on spin-isospin dependence of certain observ- 
ables. One of the limitations of such an approach is that 
only certain observables can be calculated in TDMFT. 
This limitation can be quite severe. For example, both 
the elastic and total scattering cross sections are not cal- 
culable in TDMFT. 

Witten pointed out that baryon-baryon scattering ob- 
servables are expected to possess a smooth limit only for 
momenta of order N c . For this reason we will focus on 
this regime. The regime of momenta of order N® has 
been considered elsewhere where the stress has been on 
the nucleon-nucleon potential M . 

Although we do not know the TDMFT equations for 



large N c QCD, we know its general form and can deduce 
from this the leading order spin-isospin dependence of 
certain obervables. As this form may not be immediately 
familiar, it is useful to also derive results in the context 
of TDMFT for the Skyrme model. The Skyrme model 
describes nucleons as topological solitons of a chiral field 
U and is believed to capture all of the model-independent 
large N c results of QCD. It has the virtue that the mean- 
field theory is trivially equal to the classical field theory. 
In both the Skyrme model and the Hartree treatment the 
mean-field solution for a single baryon is a "hedgehog" 
form which breaks both rotational and isorotational in- 
variance. At the mean-field level this means that there 
exist manifolds of degenerate baryon states correspond- 
ing to rotated hedgehogs. These manifolds are parame- 
terized by the three independent variables which specify 
a rotation. 

One can specify these as being the three Euler angles; 
an alternative parameterization is to specify the rotation 
in terms of a matrix A given by A = ao + a-f € SU (2), 
where ao, a%, a^ and 03 are real collective variables sat- 
isfying a constraint a\ + a ■ a = 1. The standard in- 
terpretation of these states at the full quantum level is 
that the mean-field states correspond to linear combi- 
nations of nearly degenerate states — i.e. states whose 
masses are split by 0(1/N C ). The quantized states 
can be viewed as collective wave functions associated 
with the variable parameterizing the rotations of the 
hedgehog; up to normalization constants, these collec- 
tive wave functions are simply the Wigner D matrices 
@: 9(A) = (J + l/2) 1 /2 tt-1 D I =^ ni (A) , where J and I 
are spin and isospin of the quantized state. 

In TDMFT treatments of baryon-baryon scattering, 
the initial conditions are two well-separated baryons 
moving towards each other along some axis, n (which 
in an experimental situation is the beam direction), off- 
set by an impact parameter b. The initial baryon states 
are (rotated) hedgehogs. The initial states can be pa- 
rameterized by two sets of rotation parameters — one as- 
sociated with the orientation of each initial hedgehog. 
When a quantity, 0, is calculated in TDMFT its value 
will depend on the relative orientation of the two ini- 
tial hedgehogs and b. Up to 1 /N c corrections quantities 
which are amenable to calculation in TDMFT can be 
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written as 0(A\, A%, b, h, E), where Ai and A2 represent 
the initial orientations of the two hedgehog baryons, n is 
the beam direction, E is the initial kinetic energy of each 
baryon (in the center-of-mass frame) and b is the impact 
parameter (the projection perpendicular to fi of the ini- 
tial displacement of the two particles from each other). 
The orientations, A\ and A%, are associated with partic- 
ular linear combinations of spin and isospin components 
of the baryons states. By varying these orientations one 
can extract the part of the result which is associated 
with the various spin-isospin configurations of the initial 
nucleon states. This allows us to deduce our principal re- 
sult: all observables which can be deduced from TDMFT 
must be invariant under simultaneous spin and isospin 
flips (i.e. rotations through ir/2 in both spin and isospin) 
acting on cither particle up to 1/N C corrections. 

To compute an obseravble from TDMFT requires 
that at the QCD level the observable can be ex- 
pressed directly in terms of an expectation value of some 
Heisenberg-picture operator (which we will denote as O) 
in an initial plane wave state of two nuclcons. All sen- 
sitivity to the future time evolution is contained in the 
operator. First one argues that in the p ~ N c regime 
the expectation value in these plane-wave states can be 
expressed in terms of integrals of expectation values of 
initial states consisting of well localized wave packets of 
two (rotated) hedgehog baryons heading towards each 
other and off-set by an impact parameter b. The inte- 
grations are over b and the parameters specifying the 
initial rotations of the hedgehogs. 

The integrations over b is quite standard. For p ~ N c , 
the characteristic momentum is much larger than the 
inverse size of the region of interaction. Only a small part 
of the interaction region is probed coherently ||. The 
variable b is essentially classical and the various types of 
cross sections we will use as observables can be written 
as an integral over the impact parameter which specifies 
the section of the interaction region being probed. 

The integration over the hedgehog orientations follows 
the strategy of Ref . . The mean- field theory, which is 
essentially classical in nature, accurately describes all of 
the variables except the collective ones — i.e. those de- 
grees of freedom which at the classical level can be ex- 
cited with no energetic cost. These variables must be 
quantized in order to project onto physical states with 
the correct quantum numbers. The variables associated 
with the collective isorotations of the widely-separated 
hedgehogs are collective variables. The way the quanti- 
zation is realized is simple: one expresses the observable 
at the mean-field level in terms of the A's and their time 
derivatives A. The A\^ and A\^. are dynamical variables 
which can be quantized. 

To proceed further, one must analyze the time scale 
of the interaction. The interaction takes place in a time 
of order N®: the size of the baryons is order N® and in 
this regime the velocity is also N°. Thus it takes a time 
of order N° for the two baryons to substantially overlap 



and interact. This in turn, implies that the observable 
as a function A1.2 and A\^ should be independent of N c . 
Once the observable is written as a function of Ai^ and 
A\ t 2, the next step is to make a Legendre transforma- 
tion to express the observable as a function of the A's 
and their conjugate momenta. The conjugate momenta 
can be expressed in terms of the isospin and the intrinsic 
isospin in the body-fixed frame which in hedgehog mod- 
els is the spin. However, in the process of making this 
Legendre transformation, any variables A\^ get mapped 
into h,2/1 and S\$/I. Note, that I, the moment of 
inertia, is 0(N C ) as can be seen in explicit model cal- 
culations, H, and from the model independent analysis 
H . Since in our initial states the isospin and spin are of 
order N®, it is clear that if one does a 1/N C expansion 
of the Legendre transformed operator, the leading order 
terms can depend on the A's but not on the spins or 
isospins. Thus at leading order in the 1/N C expansion 
we can neglect the spin and isospins and treat the op- 
erator as a diagonal operator in A\^- A simple way to 
understand this is that because of the large moment of 
inertia the characteristic rotational period of the hedge- 
hogs is N c and, in the large N c limit, is irrelevant during 
reactions taking place during time scales of order unity. 

The operators obtained from TDMFT are of the form 
0(Ai,A2, b, n, E). The variables A1.2, E and b are fixed 
by initial conditions in the TDMFT; b will be integrated 
over in the standard manner. The observables which we 
focus are expectation values in states with fixed initial 
nucleon spins and isospins. In experiments this is ob- 
tained by creating the given initial state multiple times 
and averaging over measurements. To compute nucleon 
matrix elements of operators which are functions of Ai t 2, 
one integrates the operator weighted by the appropriate 
collective wave functions ty(A): the expectation value of 
an operator f(Ai,A%) in a two-nucleon state with spin 
projections along an arbitrary axis rrS 1 ^ and isospin 
projections to/*- 1 ' 2 ) is 
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where / dA is the SU(2) invariant measure normalized 
so that JdA= 2ir 2 . 

The observable is a cross section associated with the 
O for nucleons with fixed initial spins and isospins, 
^ (2) . It is obtained by an integration of 
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where corrections are higher order in l/N c . 

The central result of our paper is contained in Eqs. (|l|) 
and (0). The dependence on the initial spin and isospin 
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projections of the individual nucleons of the cross sec- 
tion comes about entirely from the Wigner D matri- 
ces inside the integrals of Eq. (|l|). Note, however, that 
what is relevant is the square modulus of the D ma- 
trices. There is a well-known property of the Wigner 
D matrices (D^ )* = {—l) m ~ n D J _ m _ n which implies 

1/2 1/2 

Dm, n (A) = D_ m _ n (A) . This means that up to 

l/N c corrections the observed cross section will be un- 
changed if one simultaneously flips both the spin and the 
isospin of cither particle 1 or particle 2. This strongly 
constrains the spin and isospin dependence of the cross 
section. 

To illustrate this constraint, consider some observable 
which averages over the direction of all detected parti- 
cles. Thus, apart from the spins, the only vector left in 
the problem is the beam axis n. Suppose further that 
the observable in question is time reversal, parity even 
and isoscalar. Then from general invariance properties 
one has 

a° = A + B (<? (1) ■ 5 {2) ) + C (h- «?«) (n • <r (2) ) 
+ [Aj + B 1 (*« ■ J™) +d(n- (n ■ ct (2) )) 

x(tW-t< 2 )) (3) 

where the A, B and C coefficients are functions of energy. 
The leading order amplitude must be invariant under a 
simultaneous spin and isospin flip of either of the two 
particles. This implies that Aj = 0, Bq = 0, Co = 
at leading order, since the terms they multiply change 
signs under this transformation. For a more general ob- 
servable the same scheme may be employed: one finds 
all vectors in the problem and constructs the most gen- 
eral form the amplitude can take contracting the initial 
spins into all vectors in all ways consistent with the sym- 
metry. One then imposes the large N c rule to eliminate 
all terms which change signs under a simultaneous spin 
and isospin flip of either particle. One immediate con- 
sequence of this analysis is that all observables in this 
regime which are obtained from unpolarized measure- 
ments must be isoscalars. 

There are several fundamental limitations in the use 
of TDMFT to directly compute scattering observables 
Jl0| . TDMFT is a treatment of average quantities rather 
than any particular scattering channel. Thus, one cannot 
directly compute S-matrix elements since S-matrix ele- 
ments arc defined for particular channels. Moreover, the 
form of the initial and final states are constrained by the 
form of the trial wave function which are not rich enough 
to include the correlations necessary to get the transla- 
tionally invariant time-independent plane-wave asymp- 
totic states of the full quantum theory. Thus the initial 
states are time-dependent solutions of two lumps mov- 
ing towards each other while the final states are even 
more complicated. Clearly, since TDMFT cannot com- 
pute the S-matrix it can only give sensible information 
about quantities which average over many S-matrix el- 



ements. Such observables might, for example, be as- 
sociated with various kinds of collective flow. Moreover, 
such flow must be expressible in terms of operators which 
do not encode correlations beyond those built into the 
mean-field trial states; one can follow the flows of vari- 
ous quantum numbers but not of particular correlations 
of individual hadrons. 

It is not hard to determine which experimental ob- 
servables are accessible in mean-field theory. One sim- 
ply does, or imagines doing, a TDMFT calculation. Any 
quantity which one both knows how to compute and to 
relate to an experimentally accessible quantity is fair 
game. To make this concrete we should consider the 
following experimental observable. If we scatter two nu- 
cleons at moderately high energy (E ~ N c ), one will 
generally produce a final state with a number of mesons 
along with two baryons (and at high enough energies pos- 
sibly additional baryon-anti-baryon pairs). One can do 
the following experiment — fix a detector at angle, 9 rela- 
tive to the beam axis and measure whether an outgoing 
baryon (or anti-baryon) hits the detector. By dividing 
the number of baryons striking the detector minus the 
number of anti-baryons in solid angle dQ per unit time 
by the incident flux, one obtains a semi-inclusive par- 
tial differential cross section, da /d£l. Note that this 
quantity is semi-inclusive; while we know the position 
of one baryon, it does not specify its correlations with 
the other baryon, nor does it specify correlations with 
outgoing mesons. As such this quantity averages over 
multiple final states and is not associated with any sin- 
gle S-matrix element. (In fact, with the kinematics re- 
stricted to p ~ N c , the total amplitude for producing 
anti-baryons is both exponentially suppressed in N c Q 
and is experimentally small so one can simply use the 
semi-inclusive baryon differential cross section in place 
of the difference of the baryon from the anti-baryon. At 
higher energies one would have to use the difference of 
the two.) 

It is easy to see how to compute this in TDMFT. First 
consider using the Skyrme model at the classical level. 
While the details of the Skyrme model need not repro- 
duce QCD, it is thought that the Skyrme model cap- 
tures all of the model independent large N c results from 
QCD. The model is given in terms of the dynamics of 
a chiral field U(r,t) £ SU(2) with baryons treated as 
solitons. We can fix initial conditions U(r, 0) and U(r, 0) 
corresponding to widely separated two hedgehogs with 
orientations Ai and Ai moving with a velocity of ±vn 
(with v = \/2MnE) and with their centers offset from 
the h axis by ±6/2. Solving the equations of motion 
gives U(r, t; A\,Ai, b, h, E) where A\,A%, b, h, E param- 
eterize the initial state. In the Skyrme model the baryon 
current is given by the standard topological current, 

Bu = e -^ r Tr[{U- 1 d y U){U- 1 d a U){U- 1 d p U)\ . (4) 
Using U(r, t; Ai, A2, b, n, E) gives the baryon current as 
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a function of time and space as parameterized by the ini- 
tial conditions: B fl (f,t;Ai,A2,b,h,E). Now consider a 
large sphere of radius R centered on the nominal classi- 
cal collision point of the two solitons (which we will take 
to be the origin) . The radius R can be taken to be large 
compared to all distance scales in the problem. The av- 
erage baryon number per unit solid angle flowing at fixed 
direction (specified by polar angle 9 and azimuthal angle 
4>) out of a single collision is simply given by 
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B(Rr(9,(f)),t;A u A 2 ,b,n,E). (5) 



Integration of dn B {9 , <fi) / dil over Ai, A 2 and b as in 
Eqs.(Q) and (0) yields the physical observable 
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B(Rf(9,^),t;A 1 ,A 2 ,b,h,E) (6) 



plus corrections higher order in N c . While the proceed- 
ing derivation was done in the context of the Skyrme 
model it is clear that we would get the same structure in 
any mean-field theory including a direct time-dependent 
Hartree calculation (were that technically possible) . And 
thus the form obtained is, in fact, model-independent. 
The semi- inclusive differential cross section, da B / dS7, is 
given by the generalization of Eq. (j^) . The spin averaged 
semi-inclusive differential cross section then is given by 
Aq and is isospin independent. 

From this example it is clear how to construct other 
observables accessible in TDMFT. One can construct a 
variable measuring the fraction of total energy flowing 
out at a fixed angle. However, the most commonly used 
observables to describe nucleon-nucleon scattering — the 
total cross section and the total and differential elas- 
tic cross sections — are not amenable to mean-field treat- 
ments. Elastic cross sections are not calculable since 
they are directly associated with particular S-matrix el- 
ements. Moreover, it seems impossible to imagine any 
way in mean-field theory to extract the elastic cross sec- 
tion. The total cross section includes forward scattering 
and, as is well known, semi-classical treatments fail for 
forward scattering j(| . 

Next we consider possible corrections to the spin- 
isospin structure of Eq. (||) and its generalizations. One 
difficulty with following Witten's strategy and begin- 
ning with TDMFT as giving the leading order result is 
that we have not systematically formulated the full 1/N C 
expansion. Nevertheless, it seems apparent that while 
generic 1 /N c terms can affect the values of the various 
coefficients in Eq. (|J), they are unlikely to change the 



spin-isospin structure. The spin-isospin structure came 
about when we neglected terms sensitive to the initial 
rotational velocities of the separated hedgehogs which 
in turn implied that each explicit power of / or S came 
with a factor of l/N c . If we are considering observables 
associated with operators of good parity and isospin, the 
spin or isospin can only come in as pairs: e.g. there are 
no terms proportional to jW +1^ but only to -I^ 2 \ 
This suggests that all corrections to Eq. (||) and its gen- 
eralizations will be of order 1/N 2 . 

It is not clear whether the data on a semi-inclusive dif- 
ferential cross sections is readily accessible in a form that 
can be compared to Eq. (||). Presumably such data has 
been collected at some point. Comparisons with Eq. (||) 
would be very interesting. It would also be of inter- 
est to see if data on other observables computable from 
TDMFT can be found. 
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